Oleic Acid Could Act as a Channel Blocker in the Inhibition of nAChR: Insights from Molecular Dynamics Simulations.

The activation of the muscular nicotinic acetylcholine receptor (nAChR) produces the opening of the channel, with the consequent increase in the permeability of cations, triggering an excitatory signal. Free fatty acids (FFA) are known to modulate the activity of the receptor as noncompetitive antagonists, acting at the membrane-AChR interface. We present molecular dynamics simulations of a model of nAChR in a desensitized closed state embedded in a lipid bilayer in which distinct membrane phospholipids were replaced by two different monounsaturated FFA that differ in the position of a double bond. This allowed us to detect and describe that the cis-18:1ω-9 FFA were located at the interface between the transmembrane segments of α2 and γ subunits diffused into the channel lumen with the consequent potential ability to block the channel to the passage of ions.

Down-regulation of COX-2 activity by 1α,25(OH)2D3 is VDR dependent in endothelial cells transformed by Kaposi's sarcoma-associated herpesvirus G protein-coupled receptor.

Our previous reports showed that 1α,25-dihydroxyvitamin D3 (1α,25(OH)2D3) has antiproliferative actions in endothelial cells stably expressing viral G protein-coupled receptor (vGPCR) associated with the pathogenesis of Kaposi's sarcoma. It has been reported that COX-2 enzyme, involved in the tumorigenesis of many types of cancers, is induced by vGPCR. Therefore, we investigated whether COX-2 down-regulation is part of the growth inhibitory effects of 1α,25(OH)2D3. Proliferation was measured in presence of COX-2 inhibitor Celecoxib (10-20 µM) revealing a decreased in vGPCR cell number, displaying typically apoptotic features in a dose dependent manner similarly to 1α,25(OH)2D3. In addition, the reduced cell viability observed with 20 µM Celecoxib was enhanced in presence of 1α,25(OH)2D3. Remarkably, although COX-2 mRNA and protein levels were up-regulated after 1α,25(OH)2D3 treatment, COX-2 enzymatic activity was reduced in a VDR-dependent manner. Furthermore, an interaction between COX-2 and VDR was revealed through GST pull-down and computational analysis. Additionally, high-affinity prostanoid receptors (EP3 and EP4) were found down-regulated by 1α,25(OH)2D3. Altogether, these results suggest a down-regulation of COX-2 activity and of prostanoid receptors as part of the antineoplastic mechanism of 1α,25(OH)2D3 in endothelial cells transformed by vGPCR.

Molecular mechanisms of 33-mer gliadin peptide oligomerisation.

The proteolytic resistant 33-mer gliadin peptide is an immunodominant fragment in gluten and responsible for the celiac disease and other gluten-related disorders. Meanwhile, the primary structure of the 33-mer is associated with the adaptive immune response in celiac patients, and the structural transformation of the 33-mer into protofilaments activates a primordial innate immune response in human macrophages. This means that accumulation, oligomerisation and structural transformation of the 33-mer could be the unknown first event that triggers the disease. Herein, we reveal the early stepwise mechanism of 33-mer oligomerisation by combining multiple computational simulations, tyrosine cross-linking, fluorescence spectroscopy and circular dichroism experiments. Our theoretical findings demonstrated that the partial charge distribution along the 33-mer molecule and the presence of glutamine that favours H-bonds between the oligomers are the driving forces that trigger oligomerisation. The high content of proline is critical for the formation of the flexible PPII secondary structure that led to a ß structure transition upon oligomerisation. Experimentally, we stabilised the 33-mer small oligomers by dityrosine cross-linking, detecting from dimers to higher molecular weight oligomers, which confirmed our simulations. The relevance of 33-mer oligomers as a trigger of the disease as well as its inhibition may be a novel therapeutic strategy for the treatment of gluten-related disorders.

Orthosteric and benzodiazepine cavities of the α1ß2γ2 GABAA receptor: insights from experimentally validated in silico methods.

γ-aminobutyric acid-type A (GABAA) receptors mediate fast synaptic inhibition in the central nervous system of mammals. They are modulated via several sites by numerous compounds, which include GABA, benzodiazepines, ethanol, neurosteroids and anaesthetics among others. Due to their potential as targets of novel drugs, a detailed knowledge of their structure-function relationships is needed. Here, we present the model of the α1ß2γ2 subtype GABAA receptor in the APO state and in complex with selected ligands, including agonists, antagonists and allosteric modulators. The model is based on the crystallographic structure of the human ß3 homopentamer GABAA receptor. The complexes were refined using atomistic molecular dynamics simulations. This allowed a broad description of the binding modes and the detection of important interactions in agreement with experimental information. From the best of our knowledge, this is the only model of the α1ß2γ2 GABAA receptor that represents altogether the desensitized state of the channel and comprehensively describes the interactions of ligands of the orthosteric and benzodiazepines binding sites in agreement with the available experimental data. Furthermore, it is able to explain small differences regarding the binding of a variety of chemically divergent ligands. Finally, this new model may pave the way for the design of focused experimental studies that will allow a deeper description of the receptor.

Multiscale modelization in a small virus: Mechanism of proton channeling and its role in triggering capsid disassembly.

In this work, we assess a previously advanced hypothesis that predicts the existence of ion channels in the capsid of small and non-enveloped icosahedral viruses. With this purpose we examine Triatoma Virus (TrV) as a case study. This virus has a stable capsid under highly acidic conditions but disassembles and releases the genome in alkaline environments. Our calculations range from a subtle sub-atomic proton interchange to the dismantling of a large-scale system representing several million of atoms. Our results provide structure-based explanations for the three roles played by the capsid to enable genome release. First, we observe, for the first time, the formation of a hydrophobic gate in the cavity along the five-fold axis of the wild-type virus capsid, which can be disrupted by an ion located in the pore. Second, the channel enables protons to permeate the capsid through a unidirectional Grotthuss-like mechanism, which is the most likely process through which the capsid senses pH. Finally, assuming that the proton leak promotes a charge imbalance in the interior of the capsid, we model an internal pressure that forces shell cracking using coarse-grained simulations. Although qualitatively, this last step could represent the mechanism of capsid opening that allows RNA release. All of our calculations are in agreement with current experimental data obtained using TrV and describe a cascade of events that could explain the destabilization and disassembly of similar icosahedral viruses.

Conserved charged amino acids are key determinants for fatty acid binding proteins (FABPs)-membrane interactions. A multi-methodological computational approach.

Based on the analysis of the mechanism of ligand transfer to membranes employing in vitro methods, Fatty Acid Binding Protein (FABP) family has been divided in two subgroups: collisional and diffusional FABPs. Although the collisional mechanism has been well characterized employing in vitro methods, the structural features responsible for the difference between collisional and diffusional mechanisms remain uncertain. In this work, we have identified the amino acids putatively responsible for the interaction with membranes of both, collisional and diffusional, subgroups of FABPs. Moreover, we show how specific changes in FABPs' structure could change the mechanism of interaction with membranes. We have computed protein-membrane interaction energies for members of each subgroup of the family, and performed Molecular Dynamics simulations that have shown different configurations for the initial interaction between FABPs and membranes. In order to generalize our hypothesis, we extended the electrostatic and bioinformatics analysis over FABPs of different mammalian genus. Also, our methodological approach could be used for other systems involving protein-membrane interactions.

Echinococcus granulosus antigen B: a Hydrophobic Ligand Binding Protein at the host-parasite interface.

Lipids are mainly solubilized by various families of lipid binding proteins which participate in their transport between tissues as well as cell compartments. Among these families, Hydrophobic Ligand Binding Proteins (HLBPs) deserve special consideration since they comprise intracellular and extracellular members, are able to bind a variety of fatty acids, retinoids and some sterols, and are present exclusively in cestodes. Since these parasites have lost catabolic and biosynthetic pathways for fatty acids and cholesterol, HLBPs are likely relevant for lipid uptake and transportation between parasite and host cells. Echinococcus granulosus antigen B (EgAgB) is a lipoprotein belonging to the HLBP family, which is very abundant in the larval stage of this parasite. Herein, we review the literature on EgAgB composition, structural organization and biological properties, and propose an integrated scenario in which this parasite HLBP contributes to adaptation to mammalian hosts by meeting both metabolic and immunomodulatory parasite demands.

IFABP portal region insertion during membrane interaction depends on phospholipid composition.

Intestinal fatty acid-binding protein (IFABP) is highly expressed in the intestinal epithelium and it belongs to the family of soluble lipid binding proteins. These proteins are thought to participate in most aspects of the biology of lipids, regulating its availability for specific metabolic pathways, targeting and vectorial trafficking of lipids to specific subcellular compartments. The present study is based on the ability of IFABP to interact with phospholipid membranes, and we characterized its immersion into the bilayer's hydrophobic central region occupied by the acyl-chains. We constructed a series of Trp-mutants of IFABP to selectively probe the interaction of different regions of the protein, particularly the elements forming the portal domain that is proposed to regulate the exit and entry of ligands to/from the binding cavity. We employed several fluorescent techniques based on selective quenching induced by soluble or membrane confined agents. The results indicate that the portal region of IFABP penetrates deeply into the phospholipid bilayer, especially when CL-containing vesicles are employed. The orientation of the protein and the degree of penetration were highly dependent on the lipid composition, the superficial net charge and the ionic strength of the medium. These results may be relevant to understand the mechanism of ligand transfer and the specificity responsible for the unique functions of each member of the FABP family.

Circular dichroism and electron microscopy studies in vitro of 33-mer gliadin peptide revealed secondary structure transition and supramolecular organization.

Gliadin, a protein present in wheat, rye, and barley, undergoes incomplete enzymatic degradation during digestion, producing an immunogenic 33-mer peptide, LQLQPF(PQPQLPY)3 PQPQPF. The special features of 33-mer that provoke a break in its tolerance leading to gliadin sensitivity and celiac disease remains elusive. Herein, it is reported that 33-mer gliadin peptide was not only able to fold into polyproline II secondary structure but also depending on concentration resulted in conformational transition and self-assembly under aqueous condition, pH 7.0. A 33-mer dimer is presented as one initial possible step in the self-assembling process obtained by partial electrostatics charge distribution calculation and molecular dynamics. In addition, electron microscopy experiments revealed supramolecular organization of 33-mer into colloidal nanospheres. In the presence of 1 mM sodium citrate, 1 mM sodium borate, 1 mM sodium phosphate buffer, 15 mM NaCl, the nanospheres were stabilized, whereas in water, a linear organization and formation of fibrils were observed. It is hypothesized that the self-assembling process could be the result of the combination of hydrophobic effect, intramolecular hydrogen bonding, and electrostatic complementarity due to 33-mer's high content of proline and glutamine amino acids and its calculated nonionic amphiphilic character. Although, performed in vitro, these experiments have revealed new features of the 33-mer gliadin peptide that could represent an important and unprecedented event in the early stage of 33-mer interaction with the gut mucosa prior to onset of inflammation. Moreover, these findings may open new perspectives for the understanding and treatment of gliadin intolerance disorders.

Similar structures but different mechanisms: Prediction of FABPs-membrane interaction by electrostatic calculation.

The role of fatty acid binding proteins as intracellular fatty acid transporters may require their direct interaction with membranes. In this way different mechanisms have been previously characterized through experimental studies suggesting different models for FABPs-membrane association, although the process in which the molecule adsorbs to the membrane remains to be elucidated. To estimate the importance of the electrostatic energy in the FABP-membrane interaction, we computationally modeled the interaction of different FABPs with both anionic and neutral membranes. Free Electrostatic Energy of Binding (dE), was computed using Finite Difference Poisson Boltzmann Equation (FDPB) method as implemented in APBS (Adaptive Poisson Boltzmann Solver). Based on the computational analysis, it is found that recruitment to membranes is facilitated by non-specific electrostatic interactions. Also energetic analysis can quantitatively differentiate among the mechanisms of membrane association proposed and determinate the most energetically favorable configuration for the membrane-associated states of different FABPs. This type of calculations could provide a starting point for further computational or experimental analysis.

Prediction of the most favorable configuration in the ACBP-membrane interaction based on electrostatic calculations.

Acyl-CoA binding proteins (ACBPs) are highly conserved 10 kDa cytosolic proteins that bind medium- and long-chain acyl-CoA esters. They act as intracellular carriers of acyl-CoA and play a role in acyl-CoA metabolism, gene regulation, acyl-CoA-mediated cell signaling, transport-mediated lipid synthesis, membrane trafficking and also, ACBPs were indicated as a possible inhibitor of diazepam binding to the GABA-A receptor. To estimate the importance of the non-specific electrostatic energy in the ACBP-membrane interaction, we computationally modeled the interaction of HgACBP with both anionic and neutral membranes. To compute the Free Electrostatic Energy of Binding (dE), we used the Finite Difference Poisson Boltzmann Equation (FDPB) method as implemented in APBS. In the most energetically favorable orientation, ACBP brings charged residues Lys18 and Lys50 and hydrophobic residues Met46 and Leu47 into membrane surface proximity. This conformation suggests that these four ACBP amino acids are most likely to play a leading role in the ACBP-membrane interaction and ligand intake. Thus, we propose that long range electrostatic forces are the first step in the interaction mechanism between ACBP and membranes.